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fSJ I We have performed ^^As nuclear magnetic resonance (NMR) and nuclear quadrupole res- 

onance (NQR) on single crystalline Bai_a;Ka:Fe2As2 for x — 0.27-1. ^^As nuclear quadruple 
resonance frequency (i/q) increases linearly with increasing x. The Knight shift K in normal 
state shows Pauli paramagnetic behavior with slight temperature T dependence. The value of 
^ . K increases gradually with increasing x. By contrast, nuclear spin-lattice relaxation rate 1/Ti 

^ ' in normal state has a strong T dependence, which indicates existence of large antiferomagnetic 

(AF) spin fluctuations for all x. The T dependence of l/Ti shows a gap- like behavior below 
^ • approximately 100 K for 0.6 < a; < 0.9. These behaviors are well explained by the change of 

' band structure with expansion of hole Fermi surfaces and shrink and disappearance of electron 

C/2 , Fermi surfaces at Brillouin zone (BZ) with increasing x. The anisotropy of l/Ti, represented by 

a ratio of 1/Tiai, to 1/Tic, is always larger than 1 for all x, which indicates that the stripe-type 
AF fluctuations is dominant in this system. The K in superconducting (SC) state decreases, 
which corresponds to appearance of spin-singlet superconductivity. The T dependence of 1/Ti 
in SC state indicates multiple-SC-gap feature. A simple two gap model analysis shows that the 
larger superconducting gap gradually decreases with increasing x from 0.27 to 1 and smaller gap 
decreases rapidly and nearly vanishes for x > 0.6 where the electron pockets in BZ disappear. 
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• 1. Introduction ity, magnetic pcnetration-dcpth measurements 

O ■ Soon after the discovery of LaFcAsOi_:,F^ with su- in Ba(Fei_^Co^)2As2 have shown existence of 

^ : perconducting (SC) transition temperature Tc = 26 K nodes in underdoped and overdoped region. 1°) For 

O ■ by Kamihara et al. in 2008, iron based superconductors BaFe2(Asi_^P:,)2, angle resolved thermal conductivity 

^ ■ have been studied extensively all over the world. High- suggested the closed nodal loops located at the flat parts 

est superconducting transition temperature of 56 K of t^e electron Fermi surface,") while three dimensional 

was soon observed in iJFeAsO^'^) ("IIU", i? is a rare- nodal structure on hole Fermi surface was proposed 

earth element), which is a high discovered next to ^y theoretical calculation using fluctuation exchange 

the value in cuprate system. Many iron based super- (FLEX) approximation. 

_ conductors were observed one after another in vlFcAs In Bai_:rK^Fe2As2 (BKFA), many experiments re- 

a ; ("IIF' A is an alkaline element), Fci+^Se ("U"), and vealed appearance of multiple full SC gap around a- = 0.4 

yl£Fe2As2 ("122", AE is an alkaline-earth element). where has a maximum value of 38 K.^^ is) sign chang- 

Among them, .l^FeaAsa {AE = Ba, Ca, Sr) occupies nig s±-wave which is mediated by spin fluctuations well 

a singular position, since high quality and large single explams a lot of experimental results, T dependence 

crystals were grown. A parent compound BaFe2As2 has of spin-lattice relaxation rate 1/Ti can also be explained 

ThCr2Si2 type crystal structure and has an antiferro- by s±-wave.20.2i) However, it is pointed out theoretically 

magnetic (AF) and tetragonal to orthorhombic crystal t^at s±-wave is very fragile to nonmagnetic impurity,22) 

structure phase transition at 140 K. The orthorhombic ^^ile iron-pnictide superconductors are experimentally 

and AF phase was suppressed and superconductivity ap- ^bust against nonmagnetic impurity. No sign changing 

peared in hole doped Bai-^K^FesAsz,'') electron doped s++-wave which is mediated by orbital fluctuations and 

Ba(Fei_:,Co^)2As2,s' and BaFc2(Asi_^P^)2^'^) in which is robust against nonmagnetic impurity is another possi- 

isovalent P substitution for As acts as chemical pressure. ble candidate for the Cooper pairing in this system.23-26) 

Indeed, the direct apphcation of pressure for BaFe2As2 Hence, the SC pairing symmetry in optimum region is 

also induces superconductivity *'^) still an open question. However, there is consensus that 

As for SC gap structure, thermal conductiv- the SC gap structure itself is fully gapped one. 

By contrast, NQR and specific heat, thermal conduc- 
tivity, and magnetic penetration-depth measurements 
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for X = 1 (KFe2As2) revealed appearance of nodal SC 
gap.^^"^^-* Small angle neutron scattering (SANS) exper- 
iment has pointed out that stable isotropic hexagonal 
vortex lattice exists when magnetic fields applies paral- 
lel to crystal c axis, which supports nodes in horizontal 
direction. ^"^^ Recent muon spin relaxation measurements 
of KFe2As2 also supports the SANS results. How- 
ever, recent angle resolved photoemission spectroscopy 
(ARPES) and specific heat measurements indicate nodes 
in vertical direction. "^^■'^^^ The horizontal nodal structure 
with s-t-wave^^-* and nodal-line SC gap structure with d- 
wave'^^-' are theoretically proposed to this end compound. 
It is noteworthy that SC gap structures in BKFA are dif- 
ferent in optimally doped x — 0.4 and in heavily over- 
doped X ~ 1. 

Important feature of BKFA is a wide range real- 
ization of SC phase for 0.2 < a; < I.'') In BKFA, 
the band structure changes associated with hole dop- 
ing. For X = 0.4, cylindrical hole Fermi surfaces exist 
around F point and cylindrical electron surfaces around 
X points of the Brillouin zone (BZ). Good nesting condi- 
tion between disconnected hole and electron Fermi sur- 
faces enhances AF spin fluctuations and is favorable 
for superconductivity mediated by AF spin fluctuations. 
Within this picture, s±-wave superconductivity is ex- 
pected. ^^"^^^ Interestingly, the inelastic neutron scatter- 
ing revealed that there is incommensurate spin fluctua- 
tion Q = ((l±2(5)7r,(l±2(5)7r,0) with (5 = 0.16 for x = 1 at 
excitation energies above 3 eV^^^ where electron Fermi 
surfaces completely disappear and apparent good nest- 
ing condition between inter bands at Fermi level does not 
exit. On the other hand, there is another theoretical pro- 
posal that SC is mediated by orbital fluctuations, which 
leads to nearly orbital independent SC gap s+_|_.^'^"^^) 
Hence, there remains much work to clarify the SC gap 
symmetry in iron pnictide superconductors. Band struc- 
ture of BKFA changes associated with hole doping. The 
hole Fermi surfaces expand with increasing x, whereas 
electron Fermi surfaces shrink gradually and disappear 
for X > 0.6.^^"^^) In BKFA, Tc changes smoothly with 
increasing x. It is very important to clarify the change of 
SC gap structures in BKFA with interpolating optimally 
doped and overdoped regions. The relation between AF 
spin fluctuations and Tc should be also studied. Recently, 
we succeeded in synthesizing high quality large single 
crystals of BKFA for 0.31 < a; < 1, which provides us a 
unique opportunity to solve above problems. 

NMR and/or NQR is suitable for study of static and 
dynamical magnetic properties, and provides valuable in- 
formation on SC gap symmetry and the gap structures 
through Knight shift (K) and spin lattice relaxation rate 
(1/Ti) measurements. In this paper, we report a ^^As- 
NMR/NQR study of BKFA. 

2. Experimental 

Single crystals of BKFA {x = 0.27, 0.31, 0.39, 0.58, 
0.64, 0.69, 0.94, 1) were grown by the BaAs and KAs 
self-flux method. X-ray diffraction showed that the crys- 
tals had a tetragonal ThCr2Si2 type structure with no 
impurity phase. The detailed procedure of single crystal 
growth is basically the same described in ref. 39 which 



reports the single crystal growth of an end member of 
the system, KFe2As2. The ratio of Ba and K was de- 
termined by the energy dispersive X-ray spectroscopy. 
Experimental error of the evaluated values was within 
5%. The c-axis parameter of the single crystals was de- 
termined by X-ray diffraction analysis and the relation 
between the composition x and lattice parameter c is 
consistent with the former reported relation between the 
nominal composition x and lattice parameter by Rotter 
et al^^ 

We determined the Tc and superconducting vol- 
ume fraction of the samples with a commercial 
supcrconducting-quantum-interference-device magne- 
tometer. The Tc's of the samples are 38.5 K (a; = 0.27), 
36.5 K (x = 0.31), 38 K (x = 0.39), 30.5 K {x ^ 0.58), 
26.5 K (.T = 0.64), 20.5 K (x = 0.69), 4.5 K (x 0.94), 
3.5 K (x = 1) and the volume fraction for all the x's is 
approximately 100% except for x = 0.27. For x = 0.27, 
the fraction is approximately 60% because the crystal 
contains the phase-separated antiferromagnetically 
ordered phase^^^ as we report the experimental evidence 
of the phase separation below. We also determined 
the Tc in the NMR-measurement magnetic field by 
the change of inductance of NMR detecting coil with 
changing T. 

The NMR/NQR experiment on the ^^As nucleus (/ = 
3/2, 7/27r = 7.292 MHz/T) was carried out using phase- 
coherent pulsed NMR/NQR spectrometers and a super- 
conducting magnet between approximately 3 and 7 T. 
The measurement was performed using a ^He cryostat. 
The NMR spectra were measured by sweeping the ap- 
plied fields at a constant resonance frequency. Magnetic 
field was applied parallel to the crystal ab plane and the 
c axis. Field alignment was performed with the eye. The 
origin of the Knight shift A" = of the ^^As nucleus was 
determined by the ^'^As NMR measurement of GaAs.^^^ 
The NQR spectra were measured by sweeping the fre- 
quency in zero magnetic field. Ti was measured by a sat- 
uration recovery method at center of As-spectrum. We 
obtained Ti at a fixed frequency of 37.15 or 43.75 MHz 
with an external field of 5.06-5.09 or 5.9-5.96 T in the T 
range of 2-300 K, respectively. 

3. Results and discussion 

3. 1 NQR frequency 

In Fig. 1, we show the NQR frequency lyq against sub- 
stitution concentration x. The experimentally obtained 
z^Q was evaluated from the ^^As NQR spectra of BKFA 
as shown in the inset of Fig 1. Bars for each data point 
correspond to the full width at half maximum (FWHM) 
of the spectra. It is expected that the spectral broaden- 
ing was brought out because of the random distribution 
of Ba^^ and K"*" in the crystal. Here, zas parameter in- 
dicates z-coordinate of As site in the unit cell. 

The I'Q generally depends on temperature and these 
spectra were obtained at the different measurement 
temperatures. However, note that such temperature 
dependence of i^q in each compound is within the bar 
in Fig. 1. Hence, it is clear that the x dependence of 
the experimental vq is nearly linear. The principal axis 
of the electric field gradient is along the crystal c-axis 
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Fig. 1. (color online) Nuclear quadrupole resonance frequency 
i/Q against concentration x. BCFA, CKFA, and BKFA 
are abbreviation for Bai_a;Cs2,Fe2 As2, Cai_3;Kj;Fe2As2, and 
Bai_^Ka;Fe2As2, respectively. The experimental data of uq for 
X = were quoted from ref. 44. Closed circles are experimental 
results of i/q obtained from '^^As NQR measurements as shown 
in the inset. Open squares and diamonds are the results of band- 
structure calculation using the experimentally obtained z pa- 
rameter of As in ThCr2Si2 tetragonal structure, z parameter is 
2-coordinate of As site in the unit cell. Closed squares and dia- 
monds are the calculated results by optimizing the z parameter 
with GGA. The NQR spectra in the inset were obtained at 40 
(x = 0.39 and 0.58), 35 {x = 0.64), 30 {x = 0.69), 10 (x = 0.94), 
and 4.2 K (x = 1). 



since the As site has a local fourfold symmetry around 
the c-axis. Hence, asymmetry parameter 77 at As site 
is basically 77 = 0. However, we will discuss the possi- 
bility of finite rj in the following subsection § 3.2. The 
definition of vq at As site in these materials is as follows: 



Se'^qQ 



- 2hi{2i-iy 

where /i, eg, eQ represent the Planck constant, the elec- 
tric field gradient (EFG), and the nuclear quadrupole 
moment, respectively. The main contribution of change 
in VQ can be simply interpreted as a linear increase of 
effective ligand valency by substitution of K for Ba in 
addition to the change in lattice parameter. These two 
changes directly bring out the change in the EFG. 

This tendency is also confirmed by the systemat- 
ics obtained in the calculated electronic structures for 
BKFA. We performed the electronic structure calcula- 
tion using WIEN2K code,"'^-' where the full-potential lin- 
earized augmented-plane-wave method (FLAPW) with 
the generalized gradient approximation (GGA) for elec- 
tron correlations is used. We also used the virtual crys- 



tal approximation for alloying effect, where we made 
calculations for hypothetical materials Bai_a;CsxFc2As2 
and Cai_.^K2;Fe2As2 to sec the electronic structure of 
Bai_2;Kj:Fe2As2. Wc used experimental lattice parame- 
ters a and c reported in ref. 4. The obtained Fermi sur- 
faces using the experimental z parameter of As in the 
TliCr2Si2 tetragonal structure are consistent with Fermi 
surfaces reported by de Haas- van Alphen effect'^^-' and 
ARPES.i^'36,38) rpj^g Yermi surfaces obtained using the 
optimized value of the z parameter was not consistent for 
the end member of the system: electron Fermi surfaces 
do not disappear at the corner X point of the BZ even for 
X = \. This tendency is already reported by Singh with 
the same procedure using local density approximation. ^'^^ 
In this system changes in the lattice parameters, Fe- 
Fe interatomic distance, and the position of the As ions 
upon doping are closely interconnected with the elec- 
tronic structure. In all the calculated results, vq exhibits 
monotonic increase with x and has a same order of the 
value compared with the experimental . Of all the re- 
sults, it is surprising that the experimental i'q has a good 
agreement with the calculated for Cai_j;K2;Fc2As2 in 
the virtual crystal approximation since such agreement 
can rarely be seen in other strongly electron-correlated 
compounds. Such good agreement between experiments 
and theory suggests that the electronic structure calcu- 
lated assuming the experimental z parameters of As can 
explain the experimental electronic state of BKFA fairly 
well. 

3.2 NMR spectra 

In Fig. 2(a) and (c), we show the NMR spectra of 
BKFA with magnetic field parallel and perpendicular to 
the ah plane above Tc or Tn. In Fig. 2(b), we show the 
temperature dependent NMR spectra of x = 0.27. We 
used the NMR data for a; from ref. 44. Because ^^As 
nucleus has quantum number / = 3/2, the nuclear-spin 
Hamiltonian with an external magnetic field i/oxt is given 
by 

6 



n 



-jH,^tIz--fHc.fK-I+-^ [ill - I{I + 1)) , (2) 



{HoKt ■ K ■ I)/Hcxt = -{Kabsin9){Ix cos 9 - Iz sint 
+ {KcCos6){Ix sin 6 + Iz cost 



Iz' = Ix sin 9 + Iz cos 0, 



(3) 
(4) 

where K {Kab and K^) and 9 represent the Knight shift 
tensor (Knight shifts) and the angle between quantiza- 
tion axis and the crystal c axis, respectively. Here, we 
took the direction of the external magnetic field as prin- 
cipal axis of quantization axis and we also assumed the 
diagonalized K and the absence of in-plane anisotropy of 
Kab- The detailed procedure to obtain Kab and Kc from 
the NMR spectra is described elsewhere. 

The FWHM of the center peak for a; = 1 is about 48 
± 2 kHz. Hence, the analysis using the above expression 
works well for a; = 1 and it turned out that the actual 
field misalignment was within approximately 5% accu- 
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Fig. 2. (color online) The NMR spectra of BKFA with magnetic 
field (a) parallel and (c) perpendicular to the ab plane above Tc 
or Tn. The data for a; = was quoted from ref. 44. The measure- 
ment temperatures are 141, 100, 60, 50, 30, 25 and 4.2 K for x = 
0, 0.27, 0.39, 0.58, 0.64, 0.69 and 1, respectively. Solid lines de- 
note the fitted parameter obtained from the analysis described 
in the text, (b) The temperature dependent NMR spectra of 
X = 0.27. 



racy.^^-* But we used the following approximation for the 
other x's since the NMR spectra are much broader than 
1. 

6 



that of a; = 
H~-7(l 



Kab)HcxtIz 



1)) {H//ab), 
(5) 



U ~ -7(1 + K,)H,xtIz + ^ - /(/ + 1)) {H//c). 

(6) 

Even with the approximation, we cannot determine Kat 
for X = 0.58, 0.64, and 0.69 with sufficient accuracy be- 
cause all the center and satellite lines are broad and the 
position of the center line is affected by both the Kat 
and j/Q : the FWHM of the center peak for these com- 
position is approximately 4-8 times broader than that for 
X = 1. Therefore, the solid lines in Fig. 2(a) for x = 0.58, 
0.64, and 0.69 are just reference containing large fitting 
error. We can determine the Kc for x = 0.58, 0.64, and 
0.69 with reasonable fitting error since the position of the 



center line is affected only by the Kc and the line width 
of the center line was rather sharper: the FWHM of the 
center peak for these composition is approximately 2-3 
times broader than that for a; = 1. The satellite lines 
of X = 0.58, 0.64, and 0.69 seem to have small different 
structures. It probably arises from appearance of local 
asymmetries by K doping. Because does not affect the 
central transition, we assume rj ^ Q for the following 
discussion. 

The NMR spectra of a; = 0.27 are symmetric at 100 
and 160 K. However, the overall feature of the spectra 
slightly changes below about 60 K. In addition to the 
symmetric center and satellite lines, broad peak appears: 
the position of this peak is just below the center position. 
Below Tc this feature becomes obvious. This is because 
the paramagnetic region and the antiferromagnetically 
ordered region are phase-separated in this compound. 
This phenomenon was previously reported for polycrys- 
talline samples. ^"^^ However, recent /.tSR study of under- 
doped BKFA shows a microscopic coexistence of SC and 
AF phases without phase separation. "'^^ 

3.3 Knight shifts 

In Fig. 3(a) and (b), we show the temperature depen- 
dence of Knight shift K{T) of BKFA with magnetic field 
parallel and perpendicular to the ab plane. For all the 
data, the K(T) exhibits weak temperature dependence. 
This tendency is similar to the isovalent substitution sys- 
tem BaFe2As2-j.Pa:.'''^^ Although we could not obtain the 
Kab for x = 0.58, 0.64, and 0.69, systematic development 
of the Kc indicates that the Kab for x = 0.58, 0.64, and 
0.69 exists between x = 0.27 and 1 and that it also ex- 
hibits weak temperature dependence. With increasing x, 
K{T) at same temperature basically increases. This is at- 
tributable to the increase of density of states at the Fermi 
level with hole doping. Thus, static magnetic property is 
nearly temperature independent and shows the slight x 
dependence: spin fluctuation at around g = is not sig- 
nificant in this system. 

The hyperfine coupling constant Ai parallel or perpen- 
dicular to the afe-plane for a; = 1 was evaluated from the 
Knight shift at the As site parallel or perpendicular to 
the a&-plane and the magnetic susceptibility parallel or 
perpendicular to the a6-plane by assuming the following 
formula, 



A^0 



K,{T) = -^xt{T) . 



(7) 
(8) 

Here, Na represents Avogadro's number, xf represents 
T dependent spin susceptibility originating from 3d elec- 
tron of Fe. Xio represents T independent terms arising 
from orbital contribution. The evaluated hyperfine cou- 
pling constants and Ac originating from the transfer 
from 3d conduction electron to ^^As nuclear spin of Fe 
are -1-26(5) and -1-13(4) kOe/^e, respectively. These val- 
ues are quite similar to those obtained for x = O.^''^ This 
indicates that the hyperfine coupling constant nearly re- 
mains constant with hole doping in this system. 

In the inset of Fig. 3(b), we show the temperature de- 
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Fig. 3. (color online) The temperature dependence of Knight 
shift of BKFA with magnetic field (a) parallel and (b) perpen- 
dicular to the ab plane. The data for x = was quoted from 
ref. 44. In the inset, the temperature dependence of Kc below 
30 K is shown for x = 0.39, 0.64, 0.69. 



pendence of below 30 K for a; = 0.39, 0.64, 0.69. The 
Kc exhibits clear decrease below Tc{H) for both the sam- 
ples. Since the Knight shifts correspond to the spin sus- 
ceptibility, this result indicates spin-singlet superconduc- 
tivity in BKFA for x = 0.39, 0.64, 0.69. Combined with 
the results for x = 1,'*^) even in the heavily hole-doped 
region, the spin part of superconducting pairing symme- 
try is singlet as realized in other iron-pnictide supercon- 
ductors^^''^^'^^'^^) which are close to parent compensated 
metals. Note that we could not find any anomaly of Kc 
associated with multiple SC gap structures. 

3.4 Spin lattice relaxation rate 1/Ti in normal state 

In Fig. 4, we show the temperature dependence of spin- 
lattice relaxation rate divided by temperature l/TiT of 
^^As of BKFA in the normal state with magnetic field 
(a) parallel and (b) perpendicular to the ab plane. The 
nuclear magnetization recovery curve was fitted by the 
following double-exponential function as expected for the 
center line of the spectrum of the nuclear spin 7 = 3/2 
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Fig. 4. (color online) Temperature dependence of spin-lattice re- 
laxation rate inversely multiplied with temperature 1/T\T of 
''^As of BKFA with magnetic field (a) parallel and (b) perpen- 
dicular to the ab plane. The data for a; = was quoted from 
ref. 44. The inset shows the typical recovery curve for x = 0.69 
obtained at 1.45, 44, 267 K. 



of the ""^As nucleus,^") 
, m(t) 



7710 



0.1 exp 



t 



0.9 exp 



(9) 



where m(t) and ttiq are nuclear magnetizations after a 
time t and enough time from the NMR saturation pulse. 
In the inset of Fig. 4, we show the typical recovery curve 
for X = 0.69 obtained at 1.45, 44, 267 K. Clearly, the 
data are well fitted to the above ideal curve with a single 
Ti component. We tried to use this formula for x = 0.27 
below about 100 K under magnetic field parallel to the ah 
plane, but could not obtain good fitting results. This is 
due to more than two Ti components below this temper- 
ature because the center line shown in Fig. 2(b) consists 
of the paramagnetic and AF states. For this composi- 
tion parallel to the c axis we could obtain Ti down to 
Tc because the center line consists of only paramagnetic 
state. 

For all the composition, the \/TiT monotonically in- 
creases down to approximately 100 K. However, be- 
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low this temperature there is clear difference. While 
the 1/TiT continuously increases down to for x ~ 
0.27, 0.31, 0.39, and 1, it saturates at around 60-80 K, 
and decreases below this temperature for x = 0.58, 
0.64, 0.69. This gap-like behavior for x = 0.58, 0.64, 
0.69 is remembrance of that in electron-doped system 
Ba(Fei_2-Co2;)2As2.'*^-' In the system the temperature de- 
pendence of l/TiT is interpreted as thermal excitation 
originating from specific band structure owing to the dis- 
appearance of the hole Fermi surfaces by electron dop- 
ing 5i) In BKFA, Fermi surfaces change with increasing 
x: hole Fermi surfaces at around the F point in the BZ be- 
comes larger, and electron Fermi surfaces at around the 
X point shrinks significantly or disappears. This Fermi 
surface change is reported by recent ARPES measure- 
ments,'^®^ and it worsen the nesting condition in BKFA. 
Since the parent compound BaFe2As2 is a compensated 
metal, it is natural to consider that the gap-like behavior 
in this system has the same origin explained for electron- 
doped system Ba(Fei_a:Coa:)2As2. However, the detailed 
change of band structure in hole-doped system is some- 
what different from that in electron-doped system. In 
Ba(Fci_2;Coa;)2As2, hole Fermi surfaces disappear with 
electron doping, and their density of states near Fermi 
level also decrease with further doping. In contrast, in 
BKFA electron Fermi surfaces do not completely disap- 
pear with hole doping, but their density of states near 
Fermi level still exist with further doping until x =1. 
This is confirmed by ARPES measurements^^' 
and also by band structure calculation.^^' '^^' The recov- 
ery of increase of the 1/TiT at low temperatures in a; = 
1 is also indirect evidence for the specific band structure 
of hole-doped system. 

The AF spin fluctuation is strong near optimally doped 
region [x = 0.3-0.4). This is confirmed by our previ- 
ous NMR measurements using polycrystalline samples^^'' 
and also by neutron diffraction measurements.^^'' The AF 
spin fluctuation also exists in heavily overdoped region 
for X ~ 0.94 and 1. Because T independent dynamical 
susceptibility becomes larger with increasing x, relative 
l/TiT for a; = 0.94 and 1 is larger than that for x = 
0.3-0.4. From the temperature dependence of l/TiT, we 
can discuss the (f-depcndent spin fluctuation of system. 



>x'I(9,wo) 



(10) 



Here, and x'^{q,uj) are g-dependent hyperfine cou- 
pling constant and perpendicular component against the 
quantization axis of imaginary part of dynamical suscep- 
tibility, respectively. Based on the self consistent renor- 
malization theory (SCR) assuming two dimensional AF 
fluctuation, this general expression for 1/TiT can be 
written as follows. 

C 



1 



= Xa 



11 



Here, xo: and ^cw arc constant term of dynamical 
susceptibility. Curie constant, and Curie- Weiss tempera- 
ture, respectively, xo includes all of the temperature in- 
dependence terms of dynamical susceptibility. The Curie- 
Weiss temperature in this definition corresponds to the 



distance from the AF instability point. The value of xo 
strongly influences evaluation of ^cw- In the present 
case, we determined xo as a product of 1/TiT at 300 K 
and coefficient e{x). e{x) reflects convergence of TiT at 
high temperatures: typical values of e{x) are 0.8 and 0.55 
for a; = and 1, respectively. In order to discuss the x 
dependence of ^cw i we should perform the fltting for the 
same range of temperature. However, as discussed above, 
gap-like behavior was observed at around 50-100 K for 
X = 0.58, 0.64, and 0.69. Therefore, it is adequate to 
discuss spin fluctuation in this system using higher tem- 
perature data. We performed the fitting using the data 
above 200 K. In Fig. 5, we summarized x dependence of 
0CW obtained from the fitting and that of Tc- Because 
wc adopted the eq. (11) at high temperatures and the 
absolute value of ^cw has ambiguity due to experimen- 
tal error of Ti, the relative values and overall feature 
of the X dependence of ^cw are important to discuss the 
spin fluctuation in this system. The x dependence of ^cw 
indicates that the AF fluctuation becomes weaker with 
increasing x and that the spin fluctuation remains even 
at the end composition a; = 1. The tendency that the 
AF spin fluctuation becomes weaker with increasing x is 
attributable to the poorer nesting condition of Fermi sur- 
face with increasing x. This is consistent with recent re- 
ports on inelastic neutron scattering of BKFA."^^'^^^ The 
^cw obtained at higher temperatures corresponds to the 
potential AF fiuctuation at high energy excitation, which 
is responsible for superconductivity, compared with that 
obtained at lower temperatures. These results show that 
there is a correlation between attracting interaction of 
superconductivity and AF spin fiuctuation to some ex- 
tent. 

Although overall feature of temperature dependence 
of l/Ti is quite similar between l/Ti with magnetic 
field parallel to ab plane and that parallel to the c axis. 
However, l/Tiab{T) is always greater than l/Tic{T) 
at every temperature for all x's. In order to discuss 
this tendency which is related with the dynamical spin 
susceptibility xl*??*^); adopt the same procedure 
described in ref. 59. As described in eq. (10), l/TiT 
is related with the if-dependent hyperfine coupling 
constant and dynamical susceptibility. This equation 
can be derived by using fiuctuation-dissipation theorem. 
The combination of the hyperfine coupling constant 
and dynamical susceptibility is intrinsically brought out 
from hyperfine field at As site. By using hyperfine field, 
1/Ti can be rewritten as follows: 



+ CXD 



di((i^hf,x(i),i^hf,.(0)) 



= (/^o7n)' {\HmA^q)\^ + l^fhf,y(wo)lil3) 

Here, direction z corresponds to the direction of external 
field. Hh{,x and -ffhf.y are hyperfine fields perpendicular 
to the z direction. We assume that the hyperfine field 
H^f at As site is brought out from spin S originating 
from the nearest neighboring four Fe sites through the 
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hyperfine interaction tensor A. 



150 



fjAs 
^hf 



A = 



' Aa D Bi 
D Ab B2 
Bi B2 A, 



(14) 
(15) 



Here, Si and Bi are spin at the i-th Fe site and hyperfine 
coupling interaction tensor between the i-th Fe site and 
As site, respectively. Ai [i = a,b,c), Bi {i = 1,2), and 
D are diagonal component of A along the i direction, off 
diagonal component related with wave vector q = (tt,0) 
or (0,7r), and off diagonal component related with wave 
vector q — (tt, tt), respectively. ^'^^ Note that the wave 
vector q — (7r,0) or (0,7r) used here corresponds to the 
wave vector q = (tt, tt) which was experimentally found 
in inelastic neutron diffraction study. By using these 
equations and assuming that diagonal components of the 
tensor and spin are equivalent within ab plane {Aab = 
Aa = Ab, Sab = Sa = Sb), we obtaiu relation between 
the ratio R = Tic/Tiab and spin S. 



0.5 + 0.5 



AcSc 

Aab Sab 



R= < 



0.5 



0.5 



Sab 

~s: 



no correlation (16a) 



(7r,0) or (0,7r) (16b) 



: (tt, tt) (16c) 

As discussed in the subsection of Knight shifts, there is 
a tendency that Aab and Kab(T) are greater than Ac 
and Kc{T) for entire x's, respectively. Hence, we may 
assume that Sab is also greater than Sc- This leads that 
the value of eq. (16a) is less than 1 and that the value 
of eq. (16b) is greater than 1.5. In Fig. 6, we plotted 
temperature dependence of the ratio R. R is greater than 
1 for all cc's within error bar. This result indicates that 
spin fluctuation with the wave vector q = (tt, 0) or (0, tt) 
or wave vector near such wave vector, which corresponds 
to the wave vector q = (tt, tt) with the notation in the BZ 
for tetragonal H/mmm structure, is dominant in entire 
X composition of BKFA. Note that the current result 
anisotropy of 1 /Ti for optimally doped composition x ~ 
0.39 is consistent with the report for x ~ 0.32 by Li et 
al.^^^ Moreover, conclusion obtained in this subsection 
is consistent with the recent inelastic neutron diffraction 
studies. ^''^ 

3.5 1/Ti in superconducting state 

In Fig. 7, we show T dependence of 1/Ti for x = 0.39, 
0.58, 0.69 and 1. Each arrow indicates under external 
field. 1/Ti has no coherence peak just below Tc and de- 
creases rapidly. T dependence of l/Ti from Tc to O.STc 
is approximately proportional to for x — 0.39, T^'^ 
for X = 0.58, T^^^ for x = 0.69. The power-law behavior 
was also observed at lower temperature. The exponent 
a of power-law T" decreases with increasing x. l/Ti for 
X = 0.39 could not be obtained at lower temperatures, 
because the signal intensity dramatically decreases and 
Ti becomes longer in SC state. 
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Fig. 5. (color online) Phase diagram of BKFA and the x depen- 
dence of obtained Curie- Weiss temperature ^cw from the fitting 
described in the text. Red bold line is guide to the eye. 
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Fig. 6. (color online) Temperature dependence of the ratio of 
1/Tiab to 1/Tlc. 



l/Ti for X ~ 0.69 under 5.1 T varies proportional to T 
at lowest temperatures. With increasing external field, T 
linear region expands. And T linear behavior disappeared 
under 2.4 T. 1/Ti for x ~ 0.69 is consistent with the 
report of Zhang et al.,^"^^ although stcplike feature just 
below Tc was not observed. 

In order to understand the change of gap structure, 
we analyze l/Ti in SC state with two gap model. For a 
simplicity we assume that both gaps have same symme- 
try such as full gapped s± or line node. Here, we assume 
line node model for s-t-wave with line node or d-wave 
symmetries. The 1/Ti in SC state is expressed as, 

1 />oo 

^ ' {Nl{Ef+Ars{Ef}f{E){l-f{E)]dE, 



cx 



j=l,2 



where N'g{Ef , NPgiEf , f{E) are the density of states 
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Fig. 7. (color online) Temperature dependence of spin-lattice re- 
laxation rate 1/Ti of ^^As of BKFA. 

(DOS), the anomalous DOS arising from the coherence 
effect of Cooper paris, and the Fermi distribution func- 
tions, respectively, n,; presents the fraction of DOS of 
the i-th gap and ni+n2 = 1. This is the same procedure 
utilized in Refs. 20 and 21. 





Table 


I. Fitted parameters 






Gap type 


X 


2Ai(0)/Tc 


2A2(0)/Tc 


ni 


Si/Ai 


s± 


0.39 


9.6 


2.5 


0.76 


0.1 


s± 


0.58 


5.8 


0.92 


0.75 


0.1 


s± 


0.69 


4.8 


0.62 


0.65 


0.1 


s± 


1.0 


4.0 


0.54 


0.51 


0.1 


d-wave 


0.39 


14.2 


4.4 


0.80 




d-wave 


0.58 


8.1 


0.91 


0.75 




d-wave 


0.69 


5.6 


0.66 


0.64 




d-wave 


1.0 


4.2 


0.48 


0.50 





In Fig. 8, we show experimentally obtained 1/Ti and 
calculated 1/Ti. Each 1/Ti is normalized by the value 
at Tc in order to clarify x dependence of l/Ti. The ex- 
perimental data for x = 1 was quoted from Ref.^^'' using 
poly crystals, because we could not perform NMR/NQR 
measurements using single crystals owing to the drastic 
reduction of NMR signals in SC state. Black dashed lines 
are obtained with fully gapped s± model and brown solid 
lines with line node model. The best fitted parameters are 
listed in Tabic. I. Although wc neglected the impurity 
effect, the obtained gap parameters for x = 0.39 are con- 
sistent with the previous report on self-flux-grown single 
crystal for x = 0.32.^^'' Both models can well reproduce 
the observed l/Ti behaviors in SC state, although the 




c 

Fig. 8. (color online) Normalized temperature dependence of nor- 
malized spin-lattice relaxation rate 1/Ti at Tc of '^^As of BKFA. 
Black dashed lines and brown solid lines are calculated 1/Ti of 
fully gapped s± model and line node model, respectively. 

fitting at low temperatures is worse than that at higher 
temperatures. 

In Fig. 9, we show the x dependence of analyzed super- 
conducting gap parameter A divided by Tc of BKFA. The 
X dependence of 2Ai/Tc for both i's shows monotonous 
decrease with increasing a;. The smaller gap ^^ijT^ 
rapidly decreases from x = 0.39 to 0.58. This proba- 
bly originates from the poorer nesting condition which 
arises from the change of band structure by hole dop- 
ing. Smaller gap values above x = 0.58 were strongly 
suppressed. It corresponds to concentrations where the 
gap behavior is observed. The x dependence of l^ijTc is 
consistent with the recent results of ARPES of BKFA.^^) 
The larger gap 2Ai/rc is related with the inclination of 
1/Ti just below Tc. Compared with the x dependence 
of 2A2/rc, the X dependence of 2Ai/rc is more gradual 
and the magnitude decreases continuously. Therefore, we 
may conclude that there is no SC symmetry change in 
BKFA and that the nodal-line SC gap structure realized 
in X = 1 should be explained with the same SC sym- 
metry as that realized in optimally doped region x 
0.4. Wc may speculate that nodal-line structure emerges 
from a; 0.7 and that it develops gradually without the 
change of SC gap symmetry in this system. However, 
we cannot completely deny the possibility of the SC gap 
symmetry change at around x ^ 0.7 since the numerical 
calculation suggests that the SC condensation energy for 
c?-wavc nodal-line gap symmetry is slightly lower than 
that for fully-gapped s-wave gap symmetry and that Tc 
and the SC gap 2Ai/Tc change continuously with doping 
in such case.^^'"^'*) 

The gap symmetry cannot be determined solely with 
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Fig. 9. (color online) The x dependence of analyzed supercon- 
ducting gap parameter A divided by Tc of Bai_i,Ki,Fe2As2. 
Solid triangles, open triangles, solid diamonds, and open dia- 
monds denote 2Ai/Tc for full gap model, 2A2/TC for full gap 
model, 2^i/Tc for nodal-line gap model, and 2A2/TC for nodal- 
line gap model, respectively. Black dashed line shows zero of 
2Ai/Tc. 



NMR measurements, which can be clearly understood 
from the calculated results. To determine gap symmetry, 
not only NMR/NQR experiments but also other exper- 
iments such as specific heat, low temperature ARPES, 
are necessary as suggested in our previous report. Re- 
cent ARPES experiment of KFA observed vertical node 
on zone central hole Fermi surfaces. '^^^ Since BKFA has 
complicated band and gap structures, theoretical inter- 
pretations which can explain these results without incon- 
sistency are required. 

Finally, we discussed the origin of T linear behavior for 
X = 0.69 with increasing magnetic field at lowest mea- 
surement temperature. According to the obtained gap 
values, ratio of Ai to A2 is nearly 8. When we assume 
the upper critical field of SC (i?c2) is proportional to 
Tc(- 20 K), Hc2 of Ai would be 20 T or more than this 
value. Thus Hc2 of A2 becomes equal to or more than 
approximately 2.5 T. By this estimation it is natural to 
consider that the smaller gap is collapsed by the external 
magnetic field and that T linear behavior under higher 
magnetic fields at lower temperatures in SC state was 
observed. 

4. Conclusion 

In summary, we performed NMR/NQR measurement 
for BKFA. The x dependence of the experimental i^q 
is linear increase with increasing concentration x. This 
tendency is consistent with obtained vq by the electronic 
structure calculation using WINE2K code. 

T dependence of NMR spectra for x = 0.27 revealed 
that the paramagnetic region and antiferromanetically 
ordered region are phase-separated. Knight shift in nor- 
mal state exhibits weak temperature dependence for all 
X. With increasing x, Knight shift tends to increase due 
to the increase of DOS at the Fermi level with hole 



doping. Knight shift in SC state exhibits decease below 
Tc(iJ) for x = 0.39, 0.64, 0.69 and the spin part of SC 
paring symmetry is also singlet in the heavily hole-doped 
region. 

T dependence of 1/TiT indicates that AF spin fluc- 
tuation exist for all x. The gap-like temperature depen- 
dence for X = 0.58, 0.64, 0.69 is also observed. The gap- 
like behavior has probably the same origin explained for 
electron-doped system. In order to estimate the strength 
of spin fluctuation, we performed the fitting from 200 to 
300 K using 2D SCR theory. The x dependence of ^cw 
indicates that the spin fiuctuation becomes weaker with 
increasing x and that the spine fluctuation remains even 
at a; = 1. The spin fluctuation tends to become weaker 
with increasing x. This tendency also relates with the 
X dependence of Tc- The ratio experimental {l/Ti)ab to 
(1/Ti)c is greater than 1 in all T for all x. This indi- 
cates that the wave vector q = (tt, tt) with the notation 
in the BZ for tetragonal structure is dominant in entire 
X. Therefore, we may conclude that there is correlation 
between attracting interaction of superconductivity and 
AF spin fluctuation to some extent. 

T dependence of l/Ti in SC state has no coherence 
peak just below Tc and decreases rapidly. In order to un- 
derstand the change of gap structure, we calculate 1/Ti 
in SC state with two gap model with s± and nodal- 
line symmetry cases. Both symmetry model can well ex- 
plain in this system, although SC symmetry could not be 
determined by only NMR experiments. Both obtained 
larger gap Ai/fceTc and smaller gap lS.2/k^Tc rapidly 
decrease from x = 0.39 to 0.58. This probably originates 
from the poor nesting condition arising from the change 
of band structure. Furthermore, /S.2/k^Tc is strongly sup- 
pressed above x — 0.58. The x dependence of Ai/fceT'c 
is decreases continuously with increasing x. Therefore, it 
can be concluded that there is basically no SC symme- 
try change in BKFA. The nodal-line SC gap structure 
may be explained with the same SC symmetry in opti- 
mum doped region. According to the analysis for 1/Ti 
for X = 0.69 under magnetic flelds, we also conclude that 
the smaller gap is collapsed by the external field at lower 
temperatures. 
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